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We report on further SLAC measurements of the Askaryan effect: coherent radio emission from charge 
asymmetry in electromagnetic cascades. We used synthetic rock salt as the dielectric medium, with cascades 
produced by GeV bremsstrahlung photons at the Final Focus Test Beam. We extend our prior discovery mea- 
surements to a wider range of parameter space and explore the effect in a dielectric medium of great potential 
interest to large scale ultra-high energy neutrino detectors: rock salt (halite), which occurs naturally in high 
purity formations containing in many cases hundreds of km 3 of water-equivalent mass. We observed strong 
coherent pulsed radio emission over a frequency band from 0.2-15 GHz. A grid of embedded dual-polarization 
antennas was used to confirm the high degree of linear polarization and track the change of direction of the 
electric-field vector with azimuth around the shower. Coherence was observed over 4 orders of magnitude of 
shower energy. The frequency dependence of the radiation was tested over two orders of magnitude of UHF 
and microwave frequencies. We have also made the first observations of coherent transition radiation from the 
Askaryan charge excess, and the result agrees well with theoretical predictions. Based on these results we have 
performed a detailed and conservative simulation of a realistic GZK neutrino telescope array within a salt-dome, 
and we find it capable of detecting 10 or more contained events per year from even the most conservative GZK 
neutrino models. 



I. INTRODUCTION 

It is now widely understood that the universe becomes 
largely opaque to photons above about 100 TeV due to 
pair production on the cosmic infrared and microwave back- 
ground [1]. Stable charged baryons or leptons do not be- 
come magnetically rigid enough to propagate over intergalac- 
tic distances until their energies are so high that they also 
suffer significant losses through interactions with the cos- 
mic microwave background , through the process first noted 
by Greisen, and Zatsepin & Kuzmin (GZK) in the early 
1960's [2, 3]. 

Extragalactic astronomy in the 10 15 ~ 21 eV energy range at 
> 0.1 Gpc scales must therefore utilize other messengers, and 
neutrinos are the most likely contender. Moreover, the ex- 
istence of sources that can produce single particles with en- 
ergies approaching 1 Zeta-electron-volt (10 21 eV= 1 ZeV = 
60 Joules) provides a compelling driver toward development 
of ultra-high energy neutrino detectors, which can shed light 
not only on acceleration processes at the sources that pro- 
duce the primary cosmic rays, but on the acute problem of 
their propagation and attenuation in the intergalactic medium 
through the GZK process. This process itself yields neutrinos 
as secondaries of the pj2.7K interactions, a process first noted 
by Berezinsky & Zatsepin [4], and the resulting GZK neu- 
trino spectrum [5], peaking at EeV energies, distinctly reflects 
the cosmic ray source spectrum and distribution and provides 
unique diagnostics on the evolutionary history of the cosmic 
ray accelerators. Because the GZK neutrino flux arises from 
the aggregate intensity of all cosmic ray sources at all epochs, 
it is also potentially the strongest source of EeV neutrinos in- 
tegrated over the entire sky, and may in fact act as a "stan- 
dard candle" for EeV neutrinos. Also, given the ~ 100 TeV 
center-of-momentum energy of GZK neutrino interactions on 
nucleons at Earth, the possibility of high-energy physics ap- 



plications cannot be ignored. The combination of wide accep- 
tance for the flux models with a strong science motivation has 
made the detection of GZK neutrinos the focus of significant 
attention for all current high energy neutrino detectors. 

This goal compels us to consider detectors with target 
masses approaching 10 42 nucleons, a Teraton of mass, or 1000 
km 3 water equivalent at EeV energies [6]. Such a target mass 
must be physically accessible, and must be able to transmit in- 
formation about embedded neutrino cascades within its enor- 
mous volume to a suitable and cost-effective detector array. A 
promising approach toward this daunting task is to utilize the 
Askaryan effect [7, 8], a process which leads to strong, coher- 
ent radio pulses from such cascades. Several known dielectric 
media, with ice and rock salt the most promising at present, 
appear to have the necessary characteristics. Antarctic ice 
has been the focus of studies for this goal for several years 
now via the Radio Ice Cerenkov Experiment (RICE) [9], and 
the Antarctic Impulsive Transient Antenna (ANITA) long- 
duration balloon project will seek to exploit the effect in a 
novel approach [10]. A recent limit at extremely high neutrino 
energies has been reported via FORTE satellite observations 
of the Greenland ice sheet [11]. There is also another class of 
experiments searching for Cerenkov emission from UHE neu- 
trinos via neutrino interactions with the lunar regolith [12- 
15]. 

Rock salt, first suggested as a possible target medium by 
Askaryan [7], has also been the focus of recent efforts, and 
it appears to show equal promise with ice in regard to poten- 
tial for detector development, along with the likely advantage 
of greater accessibility [16, 17]. We report here on experi- 
ments to further understand aspects of the Askaryan process 
with a particular focus on rock salt as the medium for the cas- 
cades and radio pulse production. Our results now extend 
measurements of the effect over the entire range of cascade 
energies and radio frequency range of interest for GZK neu- 
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trino detection. We have also explored in greater detail the 
shower calorimetry aspects of radio Cerenkov measurements, 
and have made the first measurements of the vector change 
of polarization at different locations around the shower axis. 
This latter property in particular is unique to radio Cerenkov 
detection, and the combination of all these results gives us 
confidence that the road to Teraton neutrino detectors is still 
wide open along this approach. 



II. EXPERIMENTAL SETUP 

The experiment (SLAC T460) was performed at the Final 
Focus Test Beam (FFTB) facility at SLAC in June 2002. As 
in our first SLAC measurement [18], we used gamma-rays 
from a series of Aluminum radiators of different thicknesses 
to provide secondary bremsstrahlung gamma-rays for shower 
production in our rock salt target. The geometry of the target 
is shown in Fig. 1. It was built of 1.8 kg salt bricks obtained 
from Morton Salt Inc., with a total mass of about 4 metric 
tons. The bricks are slightly trapezoidal in all of their cross 
sections to accommodate the manufacturer's mold release, but 
have average dimensions ofl5x 10x6 cm. Such bricks are 
food-grade, and specified to be at least 99.5% pure sodium 
chloride. They are manufactured by compression molding of 
crystalline salt under ~ 450 tons per in 2 pressure. Their den- 
sity was measured to be 2.08 g cm~ 3 , which is about 3% less 
than the density of pure halite mineral. 

Measurements of the water content of small samples of a 
typical brick found only trace amounts of volatiles (including 
water), below 0.17%. Laboratory measurements of sodium 
chloride crystals at comparable purities give upper limits of 
< 10~ 4 on the loss tangent tan 8, which is related to the ab- 
sorption coefficient a, for small values of tan 8, by 

oc~ ve tanS (m ) 

c 

for frequency v, real part e of the dielectric permittivity, and 
speed of light c. The implied upper limits on the attenuation 
length L a = a -1 are hundreds of meters [19] in the UHF and 
cm wave regime, and our tests through of order 0.3-0.4 m of 
the salt bricks found no measurable attenuation. In the results 
reported here we assume negligible bulk absorption of the ra- 
dio emission. Paper labels attached to the the top of the bricks 
were found to be difficult to remove, and after confirming that 
they were non-conductive and had no effect on the radio trans- 
mission, they were left alone. 

To minimize the effects of the slight gaps left when stack- 
ing the salt blocks, pure (food grade) un-iodized table salt was 
used to fill in the interstices. A 2.5 cm thick polyethylene 
sheet was also mounted along the Cherenkov radiator surface 
to provide a smooth refracting surface for the external anten- 
nas and some improvement in index matching because of the 
intermediate RF index of refraction of polyethylene. 

An array of 7 by 3 (21 total) printed-circuit board (PCB) 
broadband dual-linear-polarization bowtie antennas was em- 
bedded in the upper portion of the target, on a rectangular grid 
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FIG. 1 : Geometry of the salt-block shower target used in the experi- 
ment. 



with spacings of 22 cm along the shower axis and 21 cm trans- 
verse to it, with one set of antennas arranged with one of their 
polarization axes aligned directly above the beamline. These 
are shown schematically in Fig. 1. The antennas had a fre- 
quency response of approximately 0.2-2 GHz, and about -12 
dB of cross-polarization rejection. 

In addition to the embedded array, we also used two ex- 
ternal antennas, a C/X-band standard gain horn usable over 
the frequency range from 5 to 9 GHz, at a fixed location on 
the side of the target; and a 1-18 GHz log -periodic dipole ar- 
ray (LPDA) antenna with a single linear polarization, which 
was also located to the side of the target but could be repo- 
sitioned. One side of the salt stack was arranged with a 10° 
angle with respect to the beamline to facilitate transmission 
of the Cerenkov radiation, which would have otherwise been 
totally internally reflected from a parallel face. The partially 
constructed salt stack is shown in Fig. 2 at the antenna layer, 
with several of the antennas evident, assembled with their 
leads projecting through a drilled salt brick. 

All antenna signals were digitized with Tektronix digital 
oscilloscopes (a TDS694C 3 GHz bandwidth, 10 Gsamples/s; 
and a CSA8000 sampling scope, 20 GHz, up to 50 Gsam- 
ples/s) using an ultra-stable microwave transition-radiation- 
based trigger from an upstream location. Runs were also 
taken with the bremsstrahlung radiators in and out of the beam 
at numerous times during the data collection to test for the 
presence of accelerator backgrounds, which were found to 
be negligible. The total shower energy was varied by chang- 
ing both the electron beam intensity and the thickness of the 
bremsstrahlung radiators. These radiators varied from 0.06% 
to 1.5% of a radiation length and thus extracted the same 
fraction of the electron beam energy in bremsstrahlung pho- 
tons. The electrons were then steered away from the salt target 
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FIG. 2: View of the salt stack construction at the antenna layer. An- 
tennas were mounted with their leads passing through drilled salt 
blocks as seen in several examples on the top of the stack. 



and dumped into the FFTB beam dump approximately 20 m 
downstream. The bremsstrahlung photon spectrum from the 
radiators used is very broad with a mean energy of several 
GeV. 

SLAC beam current monitors allowed us to reduce beam 
intensity to about 10 9 electrons per bunch using conventional 
methods. Below this level, closed loop control (eg. fine- 
tuning) of the beam direction was not possible. However, by 
employing our high signal-to-noise coherent microwave tran- 
sition radiation trigger, we were able to further monitor the 
beam intensity to levels an additional 2 orders of magnitude 
below this. With the combination of radiators available to us, 
the dynamic range of the showers obtained was thus about 4 
orders of magnitude in shower energy. The excellent stability 
of the FFTB beam assured that no retuning was necessary dur- 
ing our measurements, which were repeatable once the beam 
intensity was set back to the higher levels again. 

The showers produced in such experiments consist of the 
superposition of a large number of GeV electromagnetic cas- 
cades, with energies that sum up to EeV levels, but whose de- 
velopment is determined by the convolution of a GeV cascade 
and the exponential first-interaction distribution of the input 
photons. As such it cannot simulate the leading interactions 
of a cascade initiated by a single high energy particle, either 
in the longitudinal shower distributions of particle number, or 
in the exact particle species content. However, these differ- 
ences are of little consequence for radio production, since the 
Askaryan emission is dominated by the electron excess at en- 
ergies below ^100 MeV, near the region of shower maxi- 
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FIG. 3: Longitudinal distribution of the bremsstrahlung shower, for 
both total charge and charge excess as simulated by EGS4 for rock 
salt. 

mum. The charge excess itself is quite similar in either case, 
and the main quantitative differences come in only as the log- 
arithm of the ratio of the lead particle energy to the critical 
energy, which produces a more extended shower maximum 
region in the case of a shower initiated by a high energy parti- 
cle, and gives more total track-length for Cherenkov produc- 
tion. These differences are easily quantified, and in practice 
the GeV composite showers used here can be used to rigor- 
ously test models for the Askaryan radio emission. 



in. RESULTS 

To provide a baseline for evaluation of the results, we simu- 
lated the electromagnetic showers using the Electron Gamma 
Shower 4 (EGS4) Monte Carlo code [20]. EGS4 cannot di- 
rectly simulate the radio emission from the showers, but can 
assess both longitudinal and lateral development of the total 
charged particle and photon content of the shower, as well as 
the Askaryan charge excess. 

Figure 3 shows the longitudinal shower development for a 
salt target, including the charge excess. Fig. 4 shows a se- 
ries of profiles of the shower lateral distribution of the excess 
charge at different depths in cm, indicated on the plot. In each 
case the simulation is for 1000 28.5 GeV electrons initially, 
led through a 10% bremsstrahlung radiator to fully account 
for the secondary photon spectrum. No changes are required 
to the EGS4 code to develop the ^25% negative charge ex- 
cess seen; this again demonstrates that the Askaryan effect is 
generic to electromagnetic showers. 

The simulations of the shower transverse distribution were 
qualitatively confirmed by the appearance of a deposit of dis- 
colored material along a volume of salt within the stack during 
the experiment, the evidence of which is shown in Fig. 5. The 
size and shape of this cross section of the contaminated region 
is an indication of the size of the shower core itself, and it is 
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charge excess lateral distribution vs. shower depth 
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FIG. 4: Lateral distributions at various depths. Top: depths less than 
70 cm. Bottom: depths greater than 70 cm. 



still of order several cm in diameter even at a depth of ~ 15 
radiation lengths. 



A. Coherence & Absolute Field Strength. 

In reference [18], the coherence of the radiation was ob- 
served over a fairly limited range of shower energies, from 
about 3 x 10 17 eV to 10 19 eV composite energy. In our cur- 
rent measurements, we have extended the coherence measure- 
ments to a much larger range. These measurements were 
made using averaged results from several of the bowtie an- 
tennas around the region of shower maximum. Typical runs at 
a given total charge per bunch (which determines the compos- 
ite energy per bunch) were averages of 1000 beam shots, and 
the beam current was separately monitored for stability. Typi- 
cally we also recorded 100-1000 shot-to-shot pulse profiles to 
monitor the variance. Once corrected for minor variations in 
beam current, these were stable to the level of a few percent. 

Fig. 6 shows the results of the relative radio-frequency 
power as a function of shower energy, now covering nearly 
4 orders of magnitude in energy and 8 orders of magnitude in 




FIG. 5: The exit point of the beam in the downstream wall of the salt 
target, outlined by discoloration of the salt produced by beam and 
showing in cross section at the exit surface. Even at this location of 
order 15 radiation lengths along the shower, the shower core is still 
only several cm across (the salt blocks are 6 by 10 cm here in cross 
section). 
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FIG. 6: Observed coherence of the 0.3-1.5GHz radiation as a func- 
tion of total beam energy per pulse. The curve shows a quadratic 
relation for power as a function of shower energy. 
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FIG. 7: Top: Summed-in-phase pulse profile from all bowtie anten- 
nas. Bottom: Absolute power spectrum of the bowtie data, corrected 
for antenna angular response. 



RF power. No departures from the expected quadratic depen- 
dence of pulse power with shower energy are seen. This result 
demonstrates the remarkable dynamic range possible for de- 
tection of coherent RF Cerenkov radiation. 

Determination of the absolute measured field strength of 
the showers requires that the antenna effective area, spectral 
response, coupling efficiency, and angular response all be in- 
cluded in converting the measured voltages to field strengths 
in the detection medium. To properly account for the angular 
response, one must first determine if the conditions for the re- 
ception are near-field (Fresnel zone) or far-field (Fraunhofer 
Zone) for both the antenna and radiation source. The condi- 
tions for far-field response are approximately 



f< 



Rc 



L 2 cos 2 6 



(1) 



where R is the distance between source and antenna, L the 
largest antenna dimension, angle of incidence with respect 
to the normal (broadside) of the antenna, and / is the highest 
frequency for which far-field conditions are met. For the an- 
tenna as a receiving component, with R — 38 cm, L = 20 cm, 
and 9 = 25° at the Cerenkov angle the antenna far field con- 
ditions are satisfied for frequencies below 2 GHz, the highest 
usable frequency for the bowtie antennas. 

While the shower is in the far-field of the antenna, the con- 
verse is not the case, since the radiation develops over a region 
comparable in size to the entire shower with a length of about 
80 cm for full-width-at-half-maximum near the peak of cas- 
cade particle development. In practice this means that while 
the far-field antenna beam patterns may be used to determine 



the bowtie amplitude response with angle, they receive radia- 
tion only from approximately that portion of the shower that is 
projected onto the antenna effective area at the Cerenkov an- 
gle (about 66°). This behavior was observed in ref. [18], and 
provides the possibility of tracing out the longitudinal charge 
evolution of the shower by measuring the change in radiation 
amplitude along the shower axis, as was done in ref. [18] (see 
also Fig. 10 below). In determining the absolute field strength 
one must therefore account for the fact that the shower is ef- 
fectively resolved into these individual sections. 

Figure 7(top) shows a plot of the summed field strengths 
from all of the bowtie antennas, corrected for known cable at- 
tenuation but not for other frequency dependent effects. The 
pulse shape is typical of the impulse response for bowtie an- 
tennas. In the bottom portion of the figure, a Fourier ampli- 
tude spectrum of the field strength is shown, including correc- 
tions for the gain of the bowtie antennas. The solid histogram 
shows the power spectrum of the pulse, and the points are av- 
erages over 250 MHz sections of the spectrum, with the error 
bars based on the spectral variance within each section. It is 
evident that there are some near-nulls in the bowtie response, 
but the overall trend is for a rise at the higher frequencies, 
consistent with Cerenkov radiation. 

We also measured the coherence at higher frequencies (2.2- 
15.0 GHz) in several bandpasses, using the LPDA and the 
C/X-band horn antennas, both of which had known effective 
areas. Fig. 8 shows the measured electric field strength vs. 
total shower energy at different frequencies. Measurements 
were done with the LPDA on 4.95 and 14.5 GHz, respectively 
and with the horn antenna at 7.4 GHz. Points represent mea- 
sured values and lines are least-square fit curves |E| = AEj/, 01 , 
where E is the electric field and E s f, the shower energy. The 
fit for the exponent a gives: (X4.95 = 1 .00 ± 0.04 at 4.95 GHz, 
OC14.5 = 1.02 ±0.11 at 14.5 GHz, OC7.4 = 0.99 ±0.05 at 7.4 
GHz. On all three frequencies, the field strength is in agree- 
ment with full coherence of the radiation. 

Fig. 9 shows the absolute field strength measured in sev- 
eral frequency bands from 0.3-15.0 GHz. The plotted curve is 
based on the parametrization given in [21, 22], scaling from 
ice to synthetic rock salt (for details see Appendix B). The 
horizontal bars are not standard errors but show antenna or fil- 
ter bandwidth, and vertical bars indicate combined statistical 
and systematic uncertainties. The measured spectrum of the 
pulse electric field is in good agreement with the prediction. 



B. Polarization and charge excess tracking 

The crossed bowtie antennas employed in our measure- 
ments provided simultaneous measurements of the induced 
voltage in orthogonal linear polarizations, but because of their 
geometry, there is inherent leakage of a co-polarized signal 
into the cross-polarized receiver. With a measured 12 dB of 
cross-polarized rejection, up to 15% of the co-polarized am- 
plitude (eg. voltage) can appear in the cross-polarized chan- 
nel. To make accurate estimates of the angle of the projected 
plane of polarization of the radiation, the polarization leakage 
must be accounted for. We use a simple model for this, where 
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FIG. 8: Measured coherence of electric field strength at 4.95, 14.5 
with LPDA and 7.4 GHz with horn antennas respectively, with least- 
squares fit curve. 
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the co- and cross-polarized voltages are related by: 

V = E-[h + ah 90 ] (2) 
V 90 = E-[h 90 + ah ] (3) 

where E is the electric field vector of the radiation and ho,9o 
are the vector effective heights (which are complex in gen- 
eral) of the co- and cross-polarized antennas. Since each of 
the bowtie antennas are identical in both the co- and cross- 
polarized directions, and since the induced voltage has a sim- 
ple sinusoidal dependence on the projected plane of polariza- 
tion angle *P (assuming no net circular polarization), we can 
write 



V Q = E^cos^ + asinW) 
V 90 = ^(sin^ + acos 1 ?) 



(4) 
(5) 



where h — \ {\h(t\ + |/i9o|) has the same magnitude for both the 
co- and cross-polarized antennas. 

In our case we choose the co-polarized direction to be 
aligned with the beam axis, since the radiation is naturally 
expected to be polarized along this axis. If we assume that 
is approximately aligned with the beam axis (in practice this 
was done to a precision of about 1 — 2°, then for an antenna 
along the beam axis, we find to a good approximation that 
a = {V90/V0) I on-axis- is then determined by 



¥ = -tan 



_i / a-Vgo/Vo 
V0CV90/V0-I 



(6) 



Fig. 10 shows results of the polarization and relative am- 
plitude measurements for all 21 of the bowtie antennas. The 
top pane plots these as vectors with scaled lengths and direc- 
tions corresponding to the square root of the amplitude (used 
for clarity for the lower amplitudes) and the projected plane 
of the polarization l P. In the lower pane of the figure, we also 
show the relative amplitude of the antennas along the center 
axis as a function of the longitudinal distance along the shower 
(corrected for the Cerenkov angle projection). We also plot 
the EGS4 simulation of the shape of the charge excess along 
the shower (about 27% of the total charge). There is excel- 
lent agreement between the measured shape of the amplitude 
response and the normalized charge excess predictions from 
EGS4. 

In Figure 1 1 we show similar data for the three antennas at 
shower maximum at a depth of 50 cm, now plotting the plane 
of polarization as a function of the transverse position of the 
antenna with respect to the centerline of the beam. All of 
the antennas lay in a plane about 35 cm above the beamline. 
The solid curve shows the expected change in the angle of 
polarization for the three antennas, and the agreement is good. 



FIG. 9: Spectral dependence of the measured electric field strength 
for shower energy 1.9 x 10 18 eV. The solid curve is a parametrized 
Monte Carlo simulation scaled for synthetic rock salt. Vertical bars 
are estimated errors, largely due to systematics in the absolute RF 
calibration; horizontal bars indicate the bandwidth used. 



C. Transition radiation 

Transition radiation (TR) and Cerenkov radiation are 
closely related; this is particularly true in the radio regime, 
where both forms of radiation may coexist in partially or fully 
coherent forms, and interference between them may result. In 



7 




50 100 150 antenna offset from beam axis (cm) 

antenna longitudinal position (cm) 



FIG. 10: Top: Measured polarization vector field at the bowtie an- 
tenna locations in a plane 34 cm above the shower axis. The antenna 
feed center locations in that plane give the base of each vector, which 
is scaled in length as the square root of amplitude for clarity. The an- 
gles of the vectors match the observed polarization angles. Bottom: 
The observed amplitude of the received RF pulses along the center 
line of the antenna array. Also plotted is a curve of an EGS sim- 
ulation of these showers; both curves are normalized to unit peak 
amplitude. 



earlier work [23] we demonstrated the existence of free-space 
coherent microwave TR from relativistic electron bunches ex- 
iting an accelerator beam pipe, and confirmed the predicted 
angular dependence and polarization of the radiation. How- 
ever, there was a serious discrepancy in these earlier measure- 
ments with regard to the total emitted power, which was ob- 
served to be significantly lower than predictions. 

Transition radiation in the radio regime is of particular in- 
terest for detection of high energy particle showers, since 
the Askaryan charge excess guarantees not only coherent 
Cerenkov radiation, but also TR when there is any change in 
the dielectric medium on scales comparable to the size of the 
shower. Showers which break through a solid surface can pro- 
duce forward TR in the direction of propagation, and showers 
originating in one medium (air for example) which then in- 
tersect a solid surface will produce backward TR which prop- 
agates opposite to the original shower direction. It is thus of 
some importance to establish the emitted power from this pro- 
cess, to determine under what circumstances it may be useful 
for high energy shower detection. 

In a setup that was independent of the salt target, and done 
after the work with the salt was complete, we used the photon 
beam to produce an electron-gamma shower in a 2 cm block 
of lead followed by a 1.3 cm block of aluminum. The excess 



FIG. 1 1 : Plane of polarization vs. transverse distance of the antennas, 
which were all about 35 cm above the axis of the beam. 

electrons from the shower are thus expected to produce TR 
as they exit the aluminum. Strong RF pulses were observed 
with an X-band horn, with typical duration of 0.45 ns, con- 
sistent with a completely band-limited signal (h e ff = 2.4 cm, 
(A) = 4.2 cm, dv = 2 GHz, effective area A eff = 18.87 cm 2 ). 
An EGS simulation provided us with an estimate of the num- 
ber of excess electrons, their energy, and their angular and 
transverse distributions, which are shown in Fig. 12. In an ap- 
pendix, we provide details of the calculation of the expected 
TR power, which depends on the geometric form factors of 
the exiting bunch. The transverse and angular divergence pa- 
rameters are estimated to be 07 = 6 mm and Q, = 30°, respec- 
tively; therefore f T = 0.996 and % = 0.08. Most of the shower 
electrons are in the highest energy bin (0.975 < P < 1), for 
simplicity we use only the electrons in that bin for this calcu- 
lation, setting the average p = .9875. We assume the longi- 
tudinal form factor fi = 1 which is reasonable given the very 
tight longitudinal distribution of SLAC bunches. 
The results of the TR run are shown in Table I. 



TABLE I: TR measurements & analysis. 



e 


D 




N e 


Pl 


Pl 


E 


E 




cm 


msr 




mW 


mW 


V/m 


V/m 










(calc.) 


(meas.) 


(calc.) 


(meas.) 


10° 


135.5 


1.03 


5.3 x 10 8 


0.26 


0.50 


0.28 


0.31 



Here, the predicted power is about half the measured power, 
much closer than in our previous result [23]. The predicted 
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IV. DISCUSSION 
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FIG. 12: Top: Transverse spread in the simulated electron shower 
emerging from the TR radiator, grouped in 1 mm annular bins around 
the beam axis. Bottom: angular emittance spread in the same elec- 
tron shower, for the highest energy bin (0.975 < p < 1), in annular 
bins of solid angle, with 10° angular width. 



electric field is 90% of the measured electric field. We would 
expect it to be about 70% of the measured electric field based 
on the factor of 2 difference in the power, but the electric field 
estimate assumes a perfectly triangular envelope, which is not 
accurate. Detailed studies indicate that the peak and RMS 
voltages are related by Vrms = kV p /2\/2, where few 1.4 [24]. 



Our measurements have extended the validity of the theory 
underlying the Askaryan effect to a second dielectric mate- 
rial, with 30% higher density compared to silica sand, and a 
40% higher dielectric constant. In all cases the theory appears 
to scale as expected by the material properties, and there is no 
reason to believe this scaling would not apply for other dielec- 
tric materials as well. The measurements of transition radia- 
tion from the charge excess are the first observations of this 
effect, and lends further strength to the Askaryan hypothesis 
that a negative charge excess is responsible for the coherent 
radio emission, for if this were not the case, there is no obvi- 
ous way to produce the measured transition radiation with the 
strength observed. 

This latter effect could be important for both cosmic ray 
shower and neutrino shower detection, since there are many 
physical situations where a shower encounters a discontinuity 
of some kind that can lead to transition radiation. For exam- 
ple, ultra-high energy air showers encountering clouds will 
see a dielectric discontinuity; and giant air showers which im- 
pact the earth or the surface of the ocean will certainly produce 
strong backward transition radiation, particularly in the latter 
case where the Fresnel reflection coefficient from the surface 
of the ocean is of order unity. 

Our radio Cerenkov measurements demonstrate conclu- 
sively that coherent Cerenkov radio pulses from showers in 
the rock salt can accurately reconstruct the profile of the 
shower development as well as the total shower energy. They 
show that the process is inherently quadratic in the rise of 
power vs. shower energy, and this relationship is reliable 
over many decades of these parameters. Combining these 
results with in situ measurements of propagation character- 
istics at radio frequencies in salt domes [17] which indicate 
> 250 m attenuation length, we conclude that nothing in prin- 
ciple prevents the development of embedded antenna arrays 
within large salt structures such as salt domes for the charac- 
terization of ultra-high energy neutrino fluxes. 



A. Application to detector modeling. 

To further understand how this technique might be applied, 
we have created a Monte-Carlo simulation of a large-scale an- 
tenna array embedded within a generic salt dome, where the 
halite with the highest purity is likely to reside. We denote the 
array here as a Salt dome Shower Array (SalSA). The array 
is cubic with 1728 antenna nodes (12x12x12), with a grid 
spacing of 225 m, of the same order as conservative in situ es- 
timates [17] of the field attenuation length in rock salt, which 
we assume to be 250 m in this simulation. In fact recent labo- 
ratory measurements using a dielectric cavity [16] on samples 
of rock salt from the Hockly salt dome favor L a ~ 900 m at 
200 MHz, and similar long attenuation lengths are also fa- 
vored by ground-penetrating radar data [17]. 

Each node consists of 12 closely spaced antennas (0.75 m 
vertical offset), 6 of which are dipoles with vertical polariza- 
tion, and 6 of which are slotted-cylinder antennas which are 
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sensitive to horizontal polarization. The antennas are both as- 
sumed to have cos 2 dipolar response functions for simplicity 
here. The assumption is accurate for fat dipoles or bicones 
in the vertical polarization [27]; for our slot cylinders, the 
response in practice is somewhat more flattened than cos 2 
with a 1-2 dB front-to-back asymmetry, but we anticipate fur- 
ther development of current designs will tend to converge on 
a more uniform dipolar response. 

Multi-antenna local nodes such as this are effective in en- 
abling a lower local-trigger threshold in the presence of Gaus- 
sian thermal noise, and are very cost effective since the an- 
tenna elements themselves are inexpensive. The antennas are 
based on Numerical Electromagnetics (NEC2) models and 
actual tested prototypes, and are optimized for frequencies 
centered at about 200 MHz, with bandwidth of about 100%. 
Beam patterns for the mid-frequencies are assumed. 

A system temperature of 450K is assumed, based on about 
310K for the salt and a receiver noise temperature of 140K, 
consistent with low-noise amplifiers that are readily available 
commercially. The basic design for node electronics centers 
on the use of a switched-capacitor array (SCA) transient digi- 
tizer, which is only read out and sent to the surface when the 
local node is triggered; a prototype of the basic design, which 
uses Gigabit ethernet on one fiber per node to the surface, is 
described by G. Varner et al. [25, 26]. 

1. Hardware trigger rates. 

For triggering, an event must produce a voltage magnitude 
above 2.8o on 5 of 12 of the local antennas within a node. This 
threshold and coincidence level is chosen both to reduce the 
accidental rates and to ensure a robust signal for reconstruc- 
tion of the event. The accidental rates are addressed in the 
following discussion, and although we have not yet simulated 
the reconstruction process in detail, the five 2.8a signals yield 
a joint probability corresponding to 5.7a, and even assuming 
one of the five is excluded for any reason, the remaining four 
signals have a joint 5a probability. 

The propagation time across the 9m high node requires that 
the time window for such a coincidence be about 80 ns, corre- 
sponding to about 16 band-limited temporal modes of the an- 
tenna input noise voltages (x = 1/A/). For these conditions, 
random coincidence rates per node per coincidence window 
are given by the cumulative binomial probability density func- 
tion: 

where p the individual antenna probability of exceeding a 
given voltage threshold in the window, k the number of anten- 
nas required above threshold, and N T = \ N ant is the product 
of the number of antennas N ant and the number of indepen- 
dent temporal modes of length x within a coincidence window. 
Since antenna noise voltages obey a nearly Gaussian voltage 
distribution [27] 

p = -= / e- x /2 dx =0.00511. (8) 
V2k J2.8 



For bandwidths of Af ~ 200 MHz, x = 1 / Af ~ 5 ns, and thus 
the number of possible "cells" that can combine to form a 
node trigger is N z = 16x 12 = 192. The resulting probabil- 
ity of a single node trigger per 80 ns is P no de = 3.27 x 10~ 3 , 
yielding a node random trigger rate of about R noi ie = 65 kHz. 

Since this rate is too high for sustained transient digitizer 
throughput, an additional digital signal processor is used to 
reduce the rate by first-order causality constraints. To en- 
force causality for the local antenna node triggers, we re- 
quire that the hits along a node evolve along a proper light 
cone. Consider a matrix of 12 antennas by the 16 indepen- 
dent 5 ns temporal modes in the 80 ns coincidence window, 
giving a 192-cell array of antenna hits vs. time. Considering 
only plane waves for the moment, causality requires that hits 
evolve along the node on linear trajectories which must begin 
or end with the upper or lower antenna. There are ~ 16 up- 
going and <~ 16 downgoing such trajectories, each of which 
has P(p = 0.0511, fe= 5, AT X = 12) = 2.67 x 10~ 9 from equa- 
tion 7. Any of the 32 causal trajectories can begin in each 5 ns 
temporal mode, so the random rate for such node triggers is 

R n c ode ~ 2.67xlO- 9 ^p^20Hz 

which provides ample margin on the expected <~ 1 kHz 
throughput of the transient digitizer, and allows for some re- 
laxation of the trigger requirements to account for the ad- 
ditional possible patterns when wavefront curvature due to 
nearby showers is included. 

Once these local node events are transmitted to the surface, 
global triggers based on array-processing of event clusters are 
determined. An example of how this may be implemented is 
as follows. For each array node that produces a trigger, we 
require a minimum of 4 additional local node triggers among 
all of its nearest neighbors including all diagonals. For each of 
the 8 array corner nodes there are 7 neighbors, for each of the 
120 edge nodes there are 11, for each of the 600 face nodes 
there are 17, and for the remaining 1000 interior nodes there 
are 26. The light-crossing time in salt along the diagonal of an 
interior sub-cube of 27 nodes (consisting of a center node and 
its 26 neighbors) is 3.23 /js, or about 40 of the node trigger 
windows T no( i e = 80 ns. For a simple random coincidence, 
there are (27 nodes x 40 trigger windows) = 1080 cells in 
the trigger matrix. In this case, we use equation 7 with 

P = Pnode = Tnode Rnode - 80 ns X 20 Hz = 1 .6 X lO" 6 

along with k=5, and N=1080. 

The resulting cluster trigger probability for the 3.23/js win- 
dow is P duster = 1 .27 x 10~ 16 , giving a net cluster trigger rate 
of Rciuster = 3.9 x 10~ n Hz. For the same 5-node criterion the 
corresponding probabilities and rates for the faces, edges, and 
corners are far lower. Thus the potential contamination of true 
events by random coincidences per year for the entire array is 
bounded by 

Rarray = Pcluster X ^interior X 365 X 86400 = 1.23 Hz. 

indicating that, even prior to global event fitting constraints, 
the contamination from randoms can be made negligible. In 
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FIG. 13: Detection threshold estimates for a Saltdome Shower Ar- 
ray (SalSA). Solid & dashdot curves show received voltages as a 
function of distance for a 250 m (solid) and 900 m (dashdot) attenua- 
tion lengths, for the energies marked at left. The 1,2.8<T noise levels 
shown are for a 450 K system temperature. The 225 m line is the 
antenna node spacing value used for the Monte Carlo. 

practice, either local thresholds or cluster trigger criteria can 
be relaxed in hardware to ensure a manageable random "heart- 
beat" rate of order several Hz, and further constraints on clus- 
ter causality can be enforced in software. This will be done to 
ensure adequate real-time monitor of array instrument health. 



2. Thermal noise backgrounds. 

From above, causality requirements alone are adequate to 
eliminate random triggers from the neutrino candidate event 
sample. The additional information from polarization, ampli- 
tude gradient, and antenna waveform phase at each node will 
only further establish this conclusion. The strength of these 
constraints arises primarily from the fact that each node is not 
just a point detector, but provides significant directional infor- 
mation on its own. The important conclusion here is that the 
sensitivity of our implementation is limited only by our con- 
servative event reconstruction requirements (which determine 
the 5 antenna + 5 node coincidence levels), rather than random 
thermal noise backgrounds. 

We stress that thermal noise backgrounds in radio Cerenkov 
detectors, while in some ways analogous to the through-going 
cosmic -ray muon backgrounds seen in lower-energy optical 
Cerenkov detectors such as AMANDA, are in fact quite dis- 
tinct in practice, since they do not represent an irreducible 
background to the physics events, and cannot reproduce the 
characteristics of neutrino events. Atmospheric muon back- 



grounds in lower energy detectors are irreducible over some 
range of solid angle near or above the horizon. A more rele- 
vant comparison with thermal noise backgrounds can be made 
with downgoing muon events mis-reconstructed as upcom- 
ing neutrino-induced events in low-energy detectors. In this 
case, the frequency of mis-reconstructed events for the lower 
energy neutrino telescopes can be made negligible by choos- 
ing higher software thresholds, and the same is certainly true 
for radio detectors as well with regard to thermal noise back- 
grounds. 



3. Simulation results. 

Figure 13 shows the dependence of the received voltage, 
based on the parametrization of reference [21, 22], scaled for 
salt, as a function of distance and shower energy. The curves 
are labeled by shower energy in PeV, and both the la and 
2.8a levels for rms voltage above the baseline noise level are 
shown. The solid and dash-dot curves are for L a = 250, 900 m 
respectively. 

It is evident that the 225 m spacing we have employed in 
our Monte Carlo here will give a shower energy threshold in 
the neighborhood of 30 PeV for our baseline assumption of 
L a = 250 m, and in practice we find that other loss factors 
and the neutrino interaction Bjorken y-factor (eg., inelastic- 
ity) contribute to push the optimal sensitivity to neutrino en- 
ergies of order 100 PeV. Figure 13 indicates that, although it 
may be possible to operate a SalSA at a low enough energy 
threshold to detect events at the 6.4 PeV Glashow Resonance 
for W-vector boson production via v e + e~ — » W~, the node 
spacing required is below 100 m, implying an order of mag- 
nitude more array elements. Since, for a broad-spectrum neu- 
trino source, the integral over the resonance is relatively small 
in terms of the number of events [28] compared to the con- 
tinuum, the justification for the increased array density is not 
clear. In addition, the GZK neutrino spectrum peaks above 
^100 PeV. However, if L a 3> 250 m as indicated by recent 
results, the energy threshold will be considerably lower, of 
order 10 PeV at a spacing of 225 m, and only a modest in- 
crease in array density required to achieve good sensitivity at 
the Glashow resonance. 

The volume enclosed by the simulated array is approxi- 
mately 15.6 km 3 of salt, or about 34 km 3 water equivalent 
mass. An additional 11.4(63) km 3 (24.9(137) km 3 w.e.) is 
contained in the boundary of the array within one attenuation 
length (for L a = 250(900) m) of the outer wall, where we 
have indicated the perimeter values for the longer attenuation 
length in parentheses. These perimeter events will in general 
also yield events which are well measured when they illumi- 
nate a face of the array. The effective fiducial volume thus 
approaches 60(170) km 3 w.e. for an array that easily fits into 
the upper ~ 3 km of a typical salt dome. The solid angle ac- 
ceptance of the array includes the entire upper 2k hemisphere, 
and extends about 10 — 15° below the horizon, depending on 
neutrino energy. This implies of order 400(1 190) km 3 sr w.e. 
acceptance, with a threshold of order 10 eV or less. Clearly 
an accurate measure of the attenuation lengths is still an im- 



11 





FIG. 14: Four views of the same 10 18 eV hadronic shower, with directions in local altitude and azimuth (counter-clockwise from east) shown. 
The dots mark antenna nodes, and crosses mark triggered nodes, with the incoming neutrino track shown up to its termination at the shower 
vertex. 



portant factor in the design, but our goal here is to show that 
even conservative values forL a give compelling results. 

The simulation itself is a straightforward Monte Carlo in- 
tegral of the acceptance, using an isotropic neutrino flux, and 
integrating through a shell-density model for the earth to de- 
termine the angular distribution of neutrinos that can interact 
within the detector. We use a 1: 1:1 e : ft : X mix of neutrino fla- 
vors, and include neutral and charged current events. For the 
radio emission, the Zas-Alvarez-Muniz parametrization [21] 
is used. Full vector polarization is implemented. An noted 
above, we have not yet studied the event reconstruction pre- 
cision, but this is not essential to determining the sensitivity, 
once the appropriate instrumental threshold is established. 



In Figures 14, 15 we show 4 views each of two different 
events with shower energies of 10 18 and 10 19 eV. Here the 
dots locate antenna nodes, and those that produce triggers are 
marked with an 'x.' The incoming neutrino track is also in- 
dicated, with the shower occurring at its termination within 
the detector. The high refractive index of salt leads to a large 
Cerenkov angle of about 66° and accounts for the apparent 
"flatness" of the cone of triggered regions around the event 
vertices. These events are not selected to be typical but rather 
to illustrate the detector geometry. Events of lower energy, 
while able to be reconstructed via both event geometry and 
polarization information, are not as easily discerned by eye. 
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FIG. 15: Similar to previous, with shower energy = 10 eV. 



4. Estimates of GZK neutrino event rates. 

The Monte Carlo is evaluated at discrete energies over the 
PeV to ZeV energy range, assuming a monoenergetic neu- 
trino flux at each energy. The values of this function at dis- 
crete energies trace a smooth curve which represents a model- 
independent estimate of the sensitivity for fluxes which are 
smooth compared to the energy resolution of the instrument, 
which is of order AE/E ~ 1, limited primarily by the uncer- 
tainty in inelasticity, which has a mean of order 0.23 at these 
energies and a standard deviation about equal to its mean. 
The resulting curve gives the acceptance, in km 2 sr, which 
can be multiplied by the livetime to determine the exposure 
in km 2 sr s. The inverse of the exposure gives the model- 



independent flux sensitivity using the procedure outlined in 
ref. [29]. 

The result of this estimate for the flux sensitivity for 1 year 
of exposure is shown in Fig. 16, plotted along with various 
estimates of GZK neutrino fluxes [5, 30? , 31] for a complete 
range of possible parameters. A model which just equaled the 
sensitivity curve for one decade of energy would give about 
3 events. Also plotted is an earlier SalSA sensitivity esti- 
mate [33] with a simpler detector model and more conser- 
vative trigger but otherwise similar methodology. For GZK 
neutrino models we have included estimates from four inde- 
pendent calculations, and in several of these cases we provide 
both minimal and maximal estimates where they were made. 

For the three minimal GZK models shown the expected 
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FIG. 16: Broad-spectrum flux sensitivity for a SalSA for one year of 
exposure, estimated from a Monte Carlo simulation. A complete set 
of GZK models is shown, as well as an earlier estimate of SalSA sen- 
sitivity [33]. From the highest to lowest the legend indicates fluxes 
from Kalashev et al. (2002) [30] for the maximal level allowed by 
EGRET limits, an intermediate model with a 10 ZeV cutoff energy 
for the UHECR spectrum, and a minimal GZK model. The two mod- 
els marked "ESS01" are from Engel et al. 2001 [5] for one mini- 
mal model and one with strong source evolution. The chosen P&J96 
model [31] is for a minimal GZK flux. Three different models are 
shown from Aramo et al. 2004 [32], a maximal model subject only 
to observational constraints, an intermediate model which follows a 
pure E~ 2 spectrum as described by ref. [35], and another minimal 
GZK model. Finally, a model based on the unlikely assumption that 
the UHECR composition is pure iron is shown [34]. 

rates are about 11-14 events per year. If the observed rates 
are found to be significantly below this, the physics implica- 
tions for both GZK neutrinos and cosmic rays would be seri- 
ous as there are no currently self-consistent models which pre- 
dict such low rates. One class of models in which the UHECR 
are predominantly heavy nuclei such as iron have been consid- 
ered [34], and such models produce a somewhat lower GZK 
neutrino rate, leading to of order 5 events per year in our sim- 
ulated detector. The assumption of no light nuclei among the 
UHECR is however purely ad hoc at this stage, and the data 
do not support such an inference, but the importance of this 
study is that it demonstrates that there are no currently known 
ways to evade the production of GZK neutrinos, even for as- 
sumed cosmic ray compositions that are artificially designed 
to do so. 

GZK neutrino models in the mid-range that assume 
stronger source evolution or a higher-energy cutoff to the cos- 
mic ray source spectrum yield event SalSA rates in the range 
of 20-65 events per year; all of these models are still consis- 



tent with the Waxman-Bahcall limit for optically thin sources 
distributed cosmologically in a manner similar to AGN [35]. 
For the highest GZK neutrino models in which all parameters 
are maximized subject only to firm observational constraints 
such as the EGRET gamma-ray limits, the possible event rates 
are up to 120 per year. 

We note that the event rates for the minimal GZK neutrino 
flux models, while still in the regime of small statistics, are 
likely to be measured with no physics background. The only 
possible background that we currently know of comes from 
charm quark production in ultra-high energy cosmic rays, 
which can decay to energetic secondary muons, which can 
then shower via either bremsstrahlung or photonuclear pro- 
cesses in the detector. However, even the highest current esti- 
mates [36] for this process yield throughgoing muon rates of 
< 0.02 per year above 10 17 eV over a 2.4 km radius area cir- 
cumscribed around the detector volume. With regard to elec- 
tromagnetic interference, even just several meters of overbur- 
den of rock or soil above a saltdome is enough to attenuate 
all anthropogenic radio or electronic interference, so that the 
system sensitivity will be subject only to the absolute thermal 
noise floor. 

For a system with a design lifetime of 10 years or more, 
even minimal GZK neutrino fluxes will produce a rich data 
sample of events which would allow good precision on es- 
timates of the GZK neutrino energy spectrum and sky dis- 
tribution. In addition, the possibility of measuring sec- 
ondary showers from the energetic leptons in charged-current 
events compared to neutral current events would allow for 
full calorimetry of the neutral and charged-current channels, 
since at EeV energies, secondary lepton are very likely to 
shower at least 10% of their energy within 1-2 km. Landau- 
Pomeranchuk-Migdal effects on the development and radio 
emission pattern from electron-neutrino events [37] can po- 
tentially allow flavor identification, as can decays of tau lep- 
tons produced in tau neutrino events [38]. 

Dramatic increases in neutrino cross sections above the 
standard model, such as might be produced by the effects 
of microscopic black hole production in theories with large 
extra dimensions [39], would yield much larger event rates 
than described here; several thousand events per year are then 
possible even for minimal GZK models. Events so produced 
would also have quite distinct character compared to the dom- 
inant charge-current deep-inelastic neutrino events, since the 
decay of the black hole by Hawking radiation leads to a pure 
hadronic shower with no secondary high energy leptons to 
produce more extended showers. There would also be no 
significant LPM effects in electron neutrinos, assuming the 
black hole production cross section dominated over the stan- 
dard model cross sections. In fact, high cross sections com- 
bined with standard model GZK neutrino fluxes would serve 
to yield a detection with a proportionally smaller detector. 



V. CONCLUSIONS 

We conclude the following from our measurements and 
simulations: 
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1 . The Askaryan effect behaves as theoretically expected 
in measurements made in synthetic rock salt. Corre- 
lation of the strength of the radiation to the estimated 
charge excess is excellent. 

2. The polarization properties of the observed radio emis- 
sion are consistent with coherent Cerenkov radiation, 
and can be used to derive geometric properties of the 
shower track to good precision. 

3. The observed radiation is completely coherent over 
many decades of shower energy, and a wide range of 
radio and microwave frequencies. 

4. Closely related transition radiation from the Askaryan 
charge excess of showers exiting dielectric or metallic 
media has been measured and found to be in agreement 
with predictions. 

5. An antenna array of order 2.5 km on a side in an un- 
derground salt dome is well able to conclusively detect 
and characterize all current standard model GZK neu- 
trino fluxes on even a 1 year time scale; over 10 years of 
operation, good statistics can be obtained on the GZK 
neutrino spectrum above 10 17 eV. 
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of the horn. is the angle between the beamline and the line to 
the horn. £i and e 2 are the upstream and downstream complex 
dielectric constants (e^ = el + it'l), respectively. For these two 
experiments, the upstream medium is aluminum (e[ ^10) and 
the downstream medium is air (e 2 = (1.00035) 2 ). 

The energy is obtained by multiplying the RHS of equa- 
tion 1 by the bandwidth and solid angle. To convert to power, 
divide the energy by the typical time of the pulse (about the 
inverse bandwidth of the antenna). The power delivered to the 
load, Pl, is this power multiplied by the horn efficiency, about 
0.5 [27]. 

It is useful to look at the electric field E as well as the power. 
The Poynting flux S is given by S = PhlA e ff = E 2 /377, 
where A e ff is the effective area of the antenna. Therefore 
E = yJUllP^jAeff. We wish to look at the peak electric field, 

so we multiply £ by a factor of 2\/2, since E pea k = 2^/2Erms 
for a perfectly triangular envelope. 

The finite beam size introduces coherence corrections. The 
corrected power is given by 

P = N e (l+N e f L f T x)Po (11) 

where //,, fr, and % are the longitudinal, transverse and an- 
gular form factors, respectively. If the distribution of the elec- 
trons is Gaussian, then Shibata et al [40] give the form factors 
as: 

/ L = exp[-(7io L cose/A,) 2 ] (12) 



fr =exp[-(7i;c> 7 ' sinO/A,) 2 ] (13) 



/29 2 [ % l 2f> _e2, V 

X = [^2 J Q *[0 -*) K <y) + C 1 +x)E(y)]e Vdx) 

^ (14) 
where Gl and 07- are the longitudinal and transverse diver- 
gence parameters, Q. is the angular divergence parameter, and 
K(y) and E(y) are the complete elliptic integrals of the first 
and second kinds, where y = l^fxj (1 +x): 



Appendix A: Details of transition radiation analysis 

The forward spectrum of transition radiation is given in [23] 



as 



d 2 W T R _ ha A /£2sin 2 9cos 2 9 2 



d&dQ. 2ti 3 l-p 2 e 2 cos 2 
where ^ is given by 



(9) 



(e 2 - ei ) (1 - P 2 e 2 - p Vei - e 2 sin 2 9) 
(£1 + VeiVei -e 2 sin 2 9)(l - P\7ei -e 2 sin 2 9) 



(10) 



co is the bandwidth of the antenna (<f co = 2ndv), and D. is the 
solid angle subtended by the horn, given by Q. — A jd 2 where 
d is the distance from the horn to the source and A is the area 



f n / 2 

K{y)= (1 -y 2 sin 2 co)" 1/2 rfco (15) 
Jo 



Ely) = r /2 (l -y 2 sin 2 co) 1/2 Jco (16) 
Jo 

In the perfectly coherent limit, ft = fr = % = 1 and since 
N e ^l,P=N 2 P . 

To compare the measured power to the calculated power, 
we take the average of V 2 over the measured pulse and di- 
vide by 50i2. The measured voltage must be corrected for the 



attenuation in the circuit: V a 



= V m 



YQatten{dB)l20 



To get the measured electric field, take the peak electric 
field Ep ea k = 2V pea k/h e ff where h e /f is the effective height 
and Vp ea k is the measured peak voltage, with a factor of 2 to 
account for voltage division into a matched load. 
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Appendix B: Details of Cerenkov radiation analysis 

Zas, Halzen, and Stanev [42] provided a detailed analysis 
and simulations of coherent Cerenkov radio emission from 
high energy showers in ice. A variety of studies since then 
have refined these initial results and extended them to vari- 
ous scenarios of neutrino flavor and energy [43, 44]. For our 
analysis, we use the same approach as in our previous SLAC 
experiment [18], and adapt the parameterizations for the field 
strength produced by a shower in ice to the following form for 
salt: 

.ftUJI^/^Ajiy (17, 

where R is the distance from the charged-particle beam in 
Cerenkov cone direction. Aq and 8 are empirically determined 
coefficients, A = 2.53 x 1CT 7 and 8 = 1.44. The factor f d is 



a scaling factor for the difference of radiation length and den- 
sity between salt and ice, = 0.52, and K w 0.5 is a factor 
accounting for antenna near-field effects, with an estimated 
systematic uncertainty of ~ 30%. Here E s f, is the shower en- 
ergy and is given by E s h — E e N e t, where E e — 28.5 GeV is 
the energy of the electrons in the beam, N e is the number of 
electrons and t is the thickness of the radiator in Xq. The pa- 
rameters Vo = 2200 MHz, and Vi = 3500 MHz determine the 
decoherence behavior, and are estimated from the parameters 
for ice, scaled by the ratio of the radiation length in ice and 
salt, calculated from the empirical formula: 

X = 716.4 £^ — - [_£_] (18) 

where Z, is the atomic number and A; is the atomic weight of 
the elements which form the molecule. 
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